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Abstract : 

This work is devoted to the study of some exactly solvable quantum problems of four, five and six 
bodies moving on the line. We solve completely the corresponding stationary Schrodinger equation 
for these systems confined in an harmonic trap, and interacting pairwise, in clusters of two and three 
particles, by two-body inverse square Calogero potential. Both translationaly and non-translationaly 
invariant multi-body potentials are added. In each case, the full solutions are provided, namely the 
normalized regular eigensolutions and the eigenenergies spectrum. The irregular solutions are also 
studied. We discuss the domains of coupling constants for which these irregular solutions are square 
integrable. The case of a ” Coulomb-type” confinement is investigated only for the bound states of the 
four-body systems. 

PACS: 02.30.Hq, 03.65.-w, 03.65.Ge 
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1 Introduction 


A limited number of exactly solvable many-body systems exists, in the one dimensional space mm- 
The Calogero model constitutes a famous example, which was exhaustively studied mm- A survey of 
quantum integrable systems was done by Olshanetsky and Perelomov [5]. They classified the systems 
with respect to Lie algebras. Point interactions have also been considered, still in D = 1 dimensional 
space mm- 

The search for exactly solvable few-body systems concerns mostly the three-body case and still retains 
attention. Early works of three-body linear problems of Calogero-Marchioro-Wolfes [snaiio] have been 
followed by new extensions and cases. In a non exhaustive way, we can quote the three-body problems 
of Khare et al. Quesne [12] and Meljanac et al. [I3|. In a recent paper, we have exactly 

solved a generalization of the three-body Calogero-Marchioro-Wolfes problem with an additional non- 
translationally three-body potential |14| . 

Works on the four-body problem and beyond are much more scarce [151 dS). This is due to the difficulty 
to extend right away the Calogero model to a larger number of particles. We can quote the work of 
Haschke and Riihl m concerning the construction of exactly solvable quantum models of Calogero 
and Sutherland type with translationally invariant two-and four particles interactions. 

Recently an extension of the four-body problem of Wolfes m with non-translationally invariant 
interactions has been investigated and solved exactly [18]. The purpose of the present work is to point 
out particular cases admitting an exact solution. In the four-body problem, we consider the particles 
to be experiencing either an harmonic field or to be submitted to an attractive 1/r-type potential. 
The interaction among the particles is inspired by the one used in our previous work on the three-body 
potential [TT]. It contains two- and many-body forces. 

The paper is organized as follows. In section 2 we expose and solve the four-body problem for the 
case of both harmonic and ” Coulomb-type” confinement of the particles. Section 3 and 4 are devoted 
to the solutions of similar five- and six-body problems respectively. Conclusions are drawn in section 

5. 


2 A four-body problem 


2.1 The case with harmonic confinement. 


We consider the Hamiltonian 




i=l 


dx^i 


+ 




IJ- 


{Xi - XjY {xi -t X 2 -b X 3 - 3X4)2 Ya=1 


4 

Xi 


( 1 ) 
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Here, we use the units h = 2m = 1. The first contribution gives the energy of four independent 
particles in an harmonic field. To this part, residual interactions are added in the following way. The 
first three particles, with coordinates xi,X 2 and X 3 , interact pairwise by a two-body inverse square 
potentials, of Calogero type [ 8 ]. Their centre of mass interacts with the fourth particle, whose the 
coordinate is x^, via the translationally invariant term gl{x\ + X 2 + x^ — 8 x 4 )^. An additional non 
translationally invariant four-body potential, represented by the term g/{x\ X 2 + x'^-\- X 4 ), is added 
to the whole Hamiltonian. 


In order to solve this four-body problem, let us introduce the following coordinates transformation 

t = ^ixi+X2 + X3 + X 4 ), U=-^{xi-X 2 ), V = -^{xi + X2 - 2X3), W = -;^j={xi+X2 +X^-’ix^). 

( 2 ) 

The transformed Hamiltonian reads: 

^ 92 52 


dt'^ dw'^ dv? dv'^ 
, 9\{v? + v^f ^ g 

“T _ , „ _ T I 


I 2/^2 I 2 I 2 I 2 \ I 

+ ijj [t +u +v -t-rcj-l- 




+ V? + v'^ + w'^ 


2{u3-3uv^f 12tc2' 


(3) 


This Hamiltonian is not separable in {t, u, v, w} variables. To overcome this situation we introduce 
the following hyperspherical coordinates : 

t = rcosa, tc = rsinacos^, u = r sin a sin 0 sin y?, u = r sin a sin 0 cos y?, 

0 < r < oo, 0 <a< 7 r, O< 0 < 7 r, 0 <(/ 9 < 2tt. (4) 

The Schrodinger equation is then written: 


_ 3_5 2 2 , , J_ 

' Qj .2 r dr ^ ^ ^2 ^2 


52 d 

- 2 cot a— -b 
da^ da 


sin 2 a 


92 d 


-b 


-b 


92 


-b 


9A 


12cos2 0 sin2 0 I dip^ 2sin2(3(/9) 


T(r,a, 6 l, 99 ) = ET(r, a, 0 , v?) . 


(5) 


This Hamiltonian may be mapped to the problem of one particle in the four dimensional space with 
a non central potential of the form 


V{r,a,0,ip) 


fi{r) + 


1 

~9 

sin a 


f2{0) + 


sin 2 6 


( 6 ) 


It is then clear that the problem becomes separable in the four variables {r,a,6,(p}. To find the 
solution we factorize the wave function as follows : 


a, 6 *, V?) 


Vy/r sin a Vsin^ 




(7) 


3 





















Accordingly, equation ([5]) separates in four decoupled differential equations: 


(P 9A 

dif'^ 2sin^(3(/9) 


. (Bn-l) 


+ 


+ 


dO"^ ' sin^ 0 ' 12cos2 0 


0m,,n {0) — Cm,n^m,n{d^i 


d? Cm,n - i 
da^ sin^ a 


and 


, 2 2 , d‘ ^ n A 

- \- uj r A - - —^- — 

dP P 


Gl,ra,n{,Ol) — ,m,n (o:); 


( 8 ) 

(9) 

( 10 ) 

( 11 ) 


The potential of Eq.(l8|) has a periodicity of In the interval 0 < (/? < 27r, it possesses singularities 
at ip = k^,k = 0,1,...,5. This equation has been solved by Calogero [S]. In the vicinity of = 0 
(resp. ^,) the singularity can be treated if and only if A > —1/2, similar to the case of a centrifugal 
barrier. Otherwise the operator has several self-adjoint extensions, each of which may lead to a 
different spectrum [T9l [20] . 

The solutions of equation ([8|) on the interval [0,'7r/3], with Dirichlet conditions at the boundaries, in 
terms of orthogonal polynomials, are given by mm 

= (sin3(/?)5=t“C'i2'''“^(cos3vp), 0<ip<^, n = 0,1,2,... (12) 

O 

The Cn^ denote the Gegenbauer polynomials [2T|. The corresponding eigenvalues take the form 

=Q(n+^±a^ , n = 0,1,2,.... (13) 

with 

a = ^Vl + 2A . (14) 


The extension to the whole interval [0, 27r] is achieved following the prescription given in [8], by using 
symmetry arguments according to the statistics obeyed by the particles. Generally, only the regular 
solution, corresponding to 1/2 -|- a, is retained. However, the irregular solution, <l>)i \ which is 

distinct of for A > —1/2, is physically acceptable when the Dirichlet condition is satisfied for 
—1/2 < A < 0 (attractive potentials). If we release the Dirichlet condition, and ask only for the square 
integrability of the solution, as in [22|, then ^ can be retained for —1/2 < A < 3/2 [2]. For A = 0, 
the pairwise interaction between the particles 1,2 and 3 disappears. 


The second angular equation for the polar angle 9 can be written as 

' d^ . b'i-^ . /? 


62-1 

+ ^ _i + 

^ sin 2 9 


9 n Gm.n 

COS^ 9 


Qm,n{d) = 0 , 


(15) 
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( 16 ) 


where the auxiliary constants bn,f3 are defined respectively by 

~ ^ni bn = y/Sn = ^3n + 3a + —^ . 

and 

Since bn / 0, the Hamiltonian of EQ. (fTC|) is a self-adjoint operator provided that fi > —1/4. 


Considering the interval 0 < 0 < vr, we first note that Eo. dlbp has three singularities occurring at 
9 = k^, with A: = 0,1,2. This separates the interval ]0, 7 r[ in two equal length intervals namely 


0 

TT 

2 


< with r sin a cos 0 = rc > 0 , 

< 9 < TT with r sin a cos 9 = w < Q . 


(18) 

(19) 


It corresponds to a positive (or negative) value of rc, respectively. We remind the reader that sin(Q!) > 0 
as a g]0, 7 r[. In the interval 9 g]0, 7 r/ 2 [ (resp. 9 €] 7 r/ 2 , 7 r[) the centre of mass of particles 1,2 and 3 
is situated at the right (resp. the left) of the fourth particle. The presence of the centrifugal term in 
P/ cos(0)^ forbids the c.m coordinate of particles (1,2,3) to coincide with X 4 . 


The solutions of Ea. (|15l) in ]0, 7 r/ 2 [ with Dirichlet conditions are well known |14] . The regular solutions 
are 


0W(0) = (sin0)^-+5(cos0)"+^Pi/-")(cos20), 

0 < 0<'^, m = 0 ,l, 2 ,..., c= ^\/l-b4/3 



( 20 ) 

( 21 ) 


Here, are the Jacobi polynomials ^T]. The index + means that the tc-coordinate is positive. 

The solutions in the interval ^ < 9 < tt are 

0^(0) = (sin0)'>"+^(-cos0)'^+^P^’^)(cos20), (22) 

— < 9<7t, m = 0,l,2,.... 

The index — means the tc-coordinate to be negative. Note that in the interval ] 7 r/ 2 , 7 r[ we have 
— cos 9 > 0 and sin 0 > 0 so that the real power of these positive values are defined. 


In both cases, the eigenvalues Cm,n of Eq. ffTCI) are given by 

Cm,n = {2m bn P c + 1 )^, 771 = 0 , 1 , 2 ,.... (23) 


The third angular equation (fT 0 |) has been treated in HU. The regular eigensolutions and corresponding 
eigenvalues in the interval ] 0 , 7 r[ read, respectively. 


Gi, 


i,n(a) = (sin 2 2 .£ = 0,1,2,..., 


Di^r 


T Cm,n T n 


= 0 , 1 , 2 ,... . 
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(24) 

(25) 




Here, Cm,n = \/Cm,n (see Ea. (l23]l ). Note that the choice Cm,n > 0 implies for every value of {m,n} 
the function G£^m,nio() to vanish at the boundaries of the interval ]0,7r[. 


Finally, the reduced radial equation reads 


^ 2 2 , A* 4 

^ + ^2 



Fk,l,m,n{^) — 0 . 


(26) 


This is nothing but the usual 3-dimensional harmonic oscillator equation, (/r + Di^rn,n — 1/4)/r^ 
replacing the centrifugal barrier. The square integrable solutions are well known, putting a limit on 
the coefficient of the 1/r^ term. Note that taking + = 0 leads to several self-adjoint extensions 

differing by a phase. This fact has been discussed in HU. More details can be found in [231124j . It has 
to be noted that for attractive centrifugal barriers, ^ -|- D^^rn,n < 0, the problem of collapse appears, 
unless regularization procedures are taken into account [25l|26l[271|28]. Using the definition of 
Ea. (l25]i . we have 


/i -|- D^^ra,n = + {(■ + 2m -|- 3n -|- c -|- 3a -|- 3)^ >0, Vn > 0, Vm >0, Vi > 0. (27) 


The quantity /i -|- is minimal for n = 0, m = 0, ^ = 0 and a = c = 0 ( we recall that a > 0, see 

(HD and that c > 0, see (1211 1. It put constraint on //, which has to satisfy /r > —9. We introduce the 
auxiliary parameter K£^rn,n defined by 




lie,-: 




The solution of the radial equation (|26l) vanishing at r = 0 and for r ^ oo is 


(28) 


Fk,i,m,n{r) = 


exp 


ujr 


L \^(.,m,n) ( 2\ 


k = 0,1,2. 


(29) 


The are the generalized Laguerre polynomials 


The associated eigenenergies are given by 
Ek,e,m,n = 2u){2k -\- K£^m,n + 1); k = 0, 1, 2.... (30) 


Collecting all pieces, we conclude the physically acceptable solutions of the Schrodinger equation ([U) 
to be given by 




k 

TT 

o< , 


^y/fj,+{e+2m+3n+c+3a+3)'^ -lg-i^^(\/A‘+(^+2™-+3"-+c+3a+2)2 ) 

X (sin a)2-+3n+c+3a+2^(2">+3.+c+3.+3) 

X (sin (cos 29) 


X (sin 3<y9)“’'' 2 ^ ^ (cos 3<y9), 


0 , 1 , 2 ,..., 

1 — ei vr 


^ = 0,1,2,..., m = 0,l,2,...,, 


<9 < -—^-,ei = ±1, 

2 2 ““ 22 


a — —Vl + 2A 


n = 0,1,2,..., 



( 31 ) 
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Using the prescription of [8] and the parity properties of the Gegenbauer polynomials [2T], we can 
write the solution as a compact form, valid in the whole interval [0, 2ti\ for y? 




^l+(^+2m+1,n+c+‘^,a+^Y M+(^+2m+3n+c+3a+2)2 ) 

X (sin «)2m+3n+c+3a+2^(2™+3.+c+3.+3) 

X (sin 0)3’^+3“+l {ei cos (cos 26) 

xs5fn(sin3<y9)[^^“'^^^/^]| sin3(/?|“^2C'l'^''"^\cos3(^) . (32) 


Here, 62 = 1 for bosons and —1 for fermions in equation Ea. (f32|) . We recall that sgn{x) = x/\x\ 
denotes the sign of x / 0. For the Bose statistics, the extension (I52|) is possible provided that no 5 
distribution occurs when the second derivative of the wave function with respect to ip is applied at 
the boundaries between two adjacent sectors. For example, for (p = 7 r /3 and n = 0 , a 5 distribution 
occurs for a = 1/2 implying A = 0. It is due to the presence of | sin 3 (/ 9 | in (l32|) . As a consequence, the 
symmetrical solutions for the pure harmonic oscillator (A = /r = g = 0) are not recovered. Also a 6 
distribution occurs for c = 1/2 i.e. g = 0 and we do not recover the pure solutions g = 0- 


The normalization constants iV^ ^ „j,n can be calculated from 
r+oo n 2 /■(3-ei)7r/4 rf 

/ r'^dr / sin a da / sin0 d6 / dip '^k,e,m,n {r,a,9,ip)'i/k',i' ^m\n' (r,a,9,ip) 

JO Jo J{l-ei)n/4 Jo 

— ^k,k'^£,C'^m,m'^n,n'd^k,£,m,n • 

Use is made, here, of the orthogonality properties of Gegenbauer, Jacobi and Laguerre polynomials 
[29]. The full expression of the eigenenergies is expressed by 


Ek,£,m,n — Ek/-\-2m+3n — ^2fe + 1 + + (£ + 2m + 3n + c + 3o + 3)^ ^ , 

k = 0,1,2,..., £ = 0,1,2,... m = 0,l,2,... n = 0,l,2,..., . 


(34) 


Let us now consider the irregular solutions corresponding to 1/2 — a. We have to replace a by —a in 
all equations, from Ea. ffT^ to Ea. (fM)l . We recall that for —1/2 < A < 3/2, these irregular solutions 
are square integrable as mentioned before. For Fermi statistics, a 5 pathology occurs in Ea. ()32p for 
a = 1/2 (A = 0). Moreover, the requirement of self-adjointness of the Sturm-Liouville operator Ea. (|15p 
imposes us to ensure bn ^ 0 { i.e. to discard the case A = 0) and /3 > —1/4 {g > —3). Gonsequently, 
we consider values of A in ] — 1/2,0[U]0, 3/2[. 

Concerning the change a i—)• —a in the function &m,n{9), Eas. (l2UI22ll . the self-adjointness of the 
Sturm-Liouville operator imposes Bn ^ 0,Vn. Taking into account Ea. (ll3p . we have clearly to look 
for a < 1/2, which is equivalent to A < 0. Then we search for which value of a, square integrable 
solutions Sm,n{9) are obtained for every values of n. Since 

2 T = 3^ 4- 2 — 3a (o = —\/l -|- 2A) , (35) 
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we have 


(Vn >0) - + bn > 2 — 3a . (36) 

The function 0m,n(^)) Eas. (|2UI22p . leads to square integrable solutions for every n if a < 5/6. This 
last inequality happens for A < 8/9. As far as the term /3/ cos^ 0 is concerned, irregular solutions are 
found when c is changed to —c in Eas. (l20l21l22l) . These solutions are square integrable for c < 1 he. 
for /? < 3/4 or O' < 9. Consider the last angular equation concerning the variable a. The differential 
operator is self-adjoint provided that 

5 5 

(Vm > 0) (Vn > 0) Cm,n = 2m + 3n + -±c — 3a>- — 3a±c >0. (37) 

So that the self-adjointness is ensured when 

5 > 3V1 + 2A T ^1 + 1 . (38) 

The square integrability condition for the function Gi^ni,n, Ea. (l2H) . requires Cm,n + 1 > 0 Vm Vn, which 
yields 

7 

±c — 3a-|-—>0. (39) 

This defines an acceptable domain of solutions in X,g. 

As far as the radial equation is concerned, the constraint on g, + Di^ni,n > 0 is satisfied for /r > 0 
otherwise it reads 

|3 ± c - 3a| > ^/^ . (40) 

This condition defines a domain of acceptable values of /r depending on the values of A, g, within 
the interval A g] — 1/2,0[U]0,3/2[. Under such conditions, the radial solutions, Eq.([29]), are square 
integrable because K£,m,n > 0. 

Note that the spectrum for irregular solutions has eigenvalues lower than the ones corresponding to 
the regular solutions. They are given by 

= ^i^e+ 2 m+ 3 n = 1 fi + {i + 2m + 3n ± c - 3a + 3)'^ ^ , (41) 

k = 0,1,2,..., £ = 0,1,2,... m = 0,l,2,... n = 0,l,2,..., . 


2.2 The case of Coulomb-type interaction 

In their work on solvable three-body problems in H = 1, Khare and Badhuri m have considered 
an attractive interaction of the form -const ./~ spherical coordinates, such a term 

leads to an attractive 1/r potential. It can be used to replace the harmonic field in Hamiltonian like 
© and produces Coulomb-like radial wave functions. Following this idea we consider the Hamiltonian 
4 q2 Si 

^ ^ ^ (a;i+X2 + X3-3x4)2 

8 


7? > 0 (42) 
















The same coordinates transformations as Eq.(l2]) and Eq.([3]) are applied to equation Eq.(j32]), which 
becomes 



3 d 7] 1 

7 ^ ' T * <^2 1^2 


52 d 

^ - 2 cot a— + 
oa^ da 


sin^ a 


52 ^ d 


+ 


+ 


1 


92 


+ 


9A 


12cos2 0 sin2 0 I dif'^ 2sin2(399) 


^{r,a,e,ip) = E^{r,a,e,(p) . 


(43) 


Here, 'if{r,a,6,ip) are the eigensolutions associated to eigenenergy E. Compared to Eq.(IS|), the har¬ 
monic confinement has been replaced by a Coulomb-like field. 


The resolution of Ea. (l43|) is achieved in a similar manner as in the preceding section. The solutions of 
the angular parts are identical. Only the radial part is different. Using the same factorization of the 
wave function as the one of Eq.([7]), we obtain four decoupled differential equations. The eigensolutions 
and eigenvalues of the angular differential equations are given respectively by Eas. (ll2ll3D for the 
variable ip, by Eas. (l20l22p for the variable 9 and by Eas. (l24l25p for the variable a. 


The reduced radial equation reads 

/ 4 ~ d^£,m,n 4 

I j. J.2 


Ek,i,m,n 1 EfiiYn,ni'^} — O' 


(44) 


As in the case of the harmonic oscillator, the solutions for bound states ( negative energies) are well 
known. The eigenfunctions are given by [301 EH 


FkAm,n{r) = exp (-r) r) 4'"^’"*-"^(2f)r) 

V 


A: = 0,1,2..., r) = 


The associated eigenenergies take the form 


Ek,i,m,n — 


4772 


{2k -|- Ke,m,n + 1)2 


1 -|- 2k -|- 2K£ fYi j^ 


k = 0,1,2.... 


(45) 


(46) 


with Ki^rn,n given by Eq. 


In addition to the bound states, the equation Eg. (1441) admits scattering states corresponding to positive 
energies. They are not studied in this paper. 


3 A five-body problem with harmonic confinement. 


We consider the Hamiltonian 

5 

H 

i=l 





9 




















1 


9 


^ {xi - XjY (x4 - x^Y [2(a:i + X 2 + X 3 ) - 3(x4 + xY)Y 


(47) 


The first contribution gives the energy of the five independent particles in an harmonic field. Residual 
interactions are added in a following way. The first three particles, with coordinates xi,X 2 ,X 3 , consti¬ 
tuting the first cluster, interact pairwise by a two-body inverse square potentials, of Calogero type [ 8 ]. 
The remaining two-particles, whose coordinates are X 4 ,X 5 , constituting the second cluster, interact 
also via a two-body inverse square Calogero potential. The centre of mass of both clusters interact via 
the translationally invariant term g/ (2(xi +X 2 + X 3 ) — 3(x4 -|-xs))^. An additional five-body potential, 
non translationally invariant, represented by the term ^/(xf -|- X2 -|- X3 -|- X4 -|- X5), is added to the 
whole Hamiltonian. 


In order to solve this five-body problem, let us introduce the following coordinates transformation 


t = ^(xi -I- X2 -I- X3 -I- X4 -I- X5), U=^(xi-X2), X = ^(xi -I- X2 - 2 x 3 ) 

v5 v 2 v 6 

1 1 

^ =-^(x4-X5), rx = -^(2(xi -hX2 -bX3) - 3 (x 4-hxs)). 


(48) 


The transformed Hamiltonian reads: 

52 52 


H = 


52 


52 


92 


9t2 Qyfi Qy‘i 


+ oj\Y + u^ + v^ + z^ + wY 


2\2 




9\{u^ + vY 
2 {u^ — 2 >uv‘^) 


9 


2 2^2 + 30 t(;2 ■ 


(49) 


This Hamiltonian is not separable in {t, u, v, w, z} variables. As before, we introduce hyper spherical 
coordinates : 


t = rcosa, w = rsinacos9 

0 < r < 00 0 < a < vr 

z = r sin a sin 0 cos/3, u = r sin a sin 0 sin/3 sin </?, x = r sin a sin 0 sin/3 cos 
O<0<7r, 0</3<7r, 0<(/?< 27r. 


(50) 


The Schrodinger equation is then written: 

4: d 9 9 9 1 

r 

1 


^2 

dy 2 


I 2 2,9, 
-l-cxr -|— 


52 d 

-T—T - 3 cot a— -b 
oa 


sin 2 a 


92 ^ d 


-b 


9 


+ 


30 cos2 6 sin2 9 I d/3' 


92 d 

— cot + 


+ 


92 


-b 


9A 


sm 


{3 I 9y?2 2sin2(3(/?) 


9/3 2 cos2 /3 

4i{r,a,e,/3,(p) = E4>{r,a,9,P,(p) . 


(51) 
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This Hamiltonian may be mapped to the problem of one particle in the five dimensional space with a 
non central potential of the form 


V{r,a,e,ip) = /i(r) + 


9 • 2 

sin a 


f2i0) + 


sm 


fsiP) + - 


M^) \ 


sim /3 J 


(52) 


As in the previous section, the problem becomes separable but in the five variables {r, a, 6, /I, (p}. To 
find the solution we factorize the wave function as follows : 


a, „ ^ n p Gij^rn,nio‘) Qj,m,n{ 0 ) ^ . 

^k,£,j,m,nvi Oi-i 0-1 Pi 9 ■ ■ a r~-—a 

Sim/"'a smt^ Vsmp 

Accordingly, equation (IBTI) separates in five decoupled differential equations: 


(53) 






(54) 


d^ . (Bn-l) 


+ 


+ 


K 


sin"^ /3 2 cos^ /3 


— Cm,nHm,n{P') 1 


(55) 


and 


d^ {Cm,n-\) _ 

dO'^ sin^ 6 ^30 cos^ 9 


4 .+ 5 


^2 r). _ 1 ' 

_ 4 

da'^ sin^ a 


d^ I 2 2,/^”^ 4 

^ + wr H-^75- - 

dr^ 



(56) 


(57) 


(58) 


The regular solutions of equation (1541) on the interval [0,7r/3], with Dirichlet conditions at the bound¬ 
aries, are given by the expressions (see above section) 


<l>„((p) = (sin3(p)2+“C'i^'''“^(cos3(p), 0 

and associated to the eigenvalues 


vr 


<^<^1 n = 0,l,2,... 


(59) 


bl = Bn, = ( 3n + 3a-h = 0,1,2,.... a = ^Vl + 2A (60) 


The second angular equation for the polar angle (3 reads : 

( d? , bl - 


H- ^ 2 a o a ~ ) Hm^niP) — 0 . 


K 


dl3^ sin^ j3 2 cos^ /3 


(61) 
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We have bn ^ 0 so that the Hamiltonian of Ea. fjBlI) is a self-adjoint operator for k > —1/2. 

The regular solutions of Eg. (1611) in ]0,7r/2[ with Dirichlet conditions are 

= (sin/3)'’"+5(cos/3)‘'+5p^^»’'=)(cos2/3), (62) 

0 < m = 0,1,2,..., c = ^\/l 2 k , (63) 

in terms of the Jacobi polynomials. The solutions in the interval ^ < /3 < vr are 

= {sm/3f-+-^{-cos/3r+-2P^-’^\cos2l3), (64) 

^ < /3 < vr, m = 0,l,2,.... 

The index -|- (resp.—) means that the z-coordinate is positive, (resp. negative). Note that in the 
interval ]vr/2,7r[, we have — cos /3 > 0 and sin /3 > 0, so that the real power of these positive values are 
defined. 

In both cases, the eigenvalues Cm,n of Eg. ([STD are given by 

Cm,n — ^m,n^ Cm^n — (2?Ti -|- bn T C -|- 1), Tfl — 0, 1, 2, ... . (65) 


The third angular equation for the polar angle 6 can be written as 


'd9^ 


+ 


C — - 
^m,n 2 

sin^ 6 


j_i_n. 

^ on 9/1 

30cos"'6( 


&j,m,n{9) = 0, 


( 66 ) 


Since Cm,n / 0, (see Eas. (1601651) 1 the Hamiltonian of Eo. ([66|) is a self-adjoint operator provided that 
9 > -15/2. 

The regular solutions Eo. (l66p in ]0,7r/2[ with Dirichlet conditions read 

= (sin6()^"‘’"''“2 (cos 0)‘^’''2P^"*’"’‘^)(cOS 20), 

^ vr ,1 I ^ 

0 < ^ < 2 ’ "i = 0,1,2,..., ^=2V^'^T5 

The solutions in the interval ^ < 0 < vr are 

0^^(0) = (sin0)"™-"+^(-cos0)'^+^pj‘'™’"’'^^(cos20), (69) 

— < 9 <71, m = 0,1,2,.... 

The index -|- (resp.—) means the tc-coordinate to be positive, (resp. negative). Note that in the 
interval ]7r/2,7r[, we have — cos0 > 0 and sin0 > 0, so that the real power of these positive values are 
dehned. 


(67) 

( 68 ) 
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(70) 


In both cases, the eigenvalues of Ea. ([66|l are given by 

^j,m,n — ^j,m,n — (^7 Cm,n ~\~ d -\- 1), j — 0, 1, 2, ... 


The regular eigensolutions and corresponding eigenvalues of equation (I57p in the interval ]0,7r[ read, 
respectively, 


T T 2 ^ ^ 0, 1, 2, .. 


(71) 

(72) 


Note that the choice dj^rn,n > 0 implies for every value of {j,m,n} the function Gij^rn,n{ci) to vanish 
at the boundaries of the interval ]0,7r[. 


Finally, the reduced radial equation reads 


dr‘^ 


2„2 I 4 


+ + 


Ek,t,j,m,n 1 d^k,£,j,m,n{'^') — 0 . 


(73) 


Taking into account the definition of A£j^rn,n, Ea. (j72l) . we have 

/i + A£j^rn,n = /^ + (^ + 2j + 2m + 3n + c + d + 3a + 4)^ > 0 Vn > 0, Vm > 0, > 0 , (74) 

for every positive fi. The quantity fi + Aij^m,n is minimal for n = 0,j = 0,m = 0,£ = 0 and a = 0 
( we recall that a > 0, see pUjl and that c > 0, see pSIl l. It put constraint on negative values of /r, 
namely —16 < /r < 0. 


We introduce the auxiliary parameter K.£j^rn,n defined by 

^t,j,m,n — T Alj^rn,n i ^£,j,m,n — T 

The solution of the radial equation ([751) is 


FkAj,m,nir) = exp j k = 0,l, 2 ... . 

in terms of the generalized Laguerre polynomials. The associated eigenenergies are given by 

2cj(2A) T T 1)? ^ 0? f; 2 — 


(76) 


(77) 
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Collecting all pieces, we conclude the physically acceptable regular solutions of the Schrodinger equa¬ 
tion ()5ip to be given in a compact form ( as explained in the previous section) by 


1 Oij ^1 Pj 


k 


^V/.+ (£+2i+2m+3n+c+d+3a+4)= >^+iWj+2m+3n+c+d+3a+3f ) 

X (sin «)2i+2r„+3n+c+d+3a+|^W+2^+3.+c+.+4) 

X (sin 0)2m+3n+3a+c+2^^^ p{‘^rn+3n+c+3a+l,d) ^0) 

X (sin /3)3"+3“+I (ei cos (cos 2/3) 

xs(7n(sin3/3)f*'^“^3)/2] | sin3(/5|“+2C'^“^2)^^Qg 


0,1,2,..., ^ = 0,1,2,..., j = 0,l,2,.., m = 0,l,2,.., n = 0,l,2,..., 

vr 1 — eiTT ^ 3 — eiTT 

2-''==^' 

2 2 ~ ~ 2 2 


0 < a < TT 


1 


1 


= ^Vl + 2A , c= ^VI + 2k d=^Jl+‘^ . 
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Here, 63 = 1 for bosons and —1 for fermions in equation Ea. (|78l) . For the Bose statistics, the extension 
(j78p is possible provided that no 5 distribution appears when the second derivative of the wave function 
is applied at the boundaries between two adjacent sectors. For cp = tt/3 and n = 0, a <5 distribution 
occurs for a = 1/2 implying A = 0, due to the term | sin 3 (^ 9 1 in (j78p . Also a 6 distribution occurs 
for d = 1/2 i.e. g = 0 at 9 = 7r/2. Consequently, the symmetrical solutions for the pure harmonic 
oscillator {X = g = g = 0) are not recovered. 


The normalization constants can be calculated from 

^- 1-00 . „ r(3-e2)vi-/4 r(3-ei)7r/4 

/ r^dr / sim a da sim 6 d9 sin /3 d/3 

Jo Jo </(l— 62 ) 77/4 </(!— 61 ) 77/4 

TT 

X / dip'^£ a, 9, P, £/ji jYi'P (^T, a, 9, P, ^k,k'^i,l'^j,j'^m,m'^n,n'^k,l,j,m,n ■('^9) 

JO 


As above, use is made of the orthogonality properties of Gegenbauer, Jacobi and Laguerre polynomials. 


The full expression of the eigenenergies is expressed by 


Ek,i,j,m,n — Ek/-\-2j+2m+3n — ^2fc + 1 + \Jg, -\- {i -\- 2j -\- 2m + 3n + c + d + 3o -|- 4)^ ^ 

k = 0,1,2,..., ^ = 0,1,2,... j = 0,1,2, m = 0,1,2,... n = 0,1,2,.... (80) 
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4 A six-body problem with harmonic confinement. 


We consider the Hamiltonian 


6/^2 \ ^ I ^ I 

^ +^iE7-T2+^2 E 7-72 


i=l 


+ 


9 


+ 


9 


(xi + X2 + X3 - X4 - X5 - Xe)^ Ta= 1 xj 


(81) 


The first contribution gives the energy of the six independent particles in an harmonic field. The first 
three particles, whose coordinates are xi,X 2 ,X 3 , constitute the first cluster and the remaining three 
particles, whose coordinates are X 4 ,X 5 ,X 6 , constitute the second cluster. In each cluster, the particles 
interact pairwise by a two-body inverse square potentials. The center of mass of both clusters interact 
via the six-body translationally invariant inverse square potential, namely g/ (X1-I-X2-I-X3—X4—X5— xe)^. 
A non-translationaly invariant six-body potential with coupling constant g is added. 

In order to solve this six-body problem, let us introduce the following coordinates transformation 
t = _l_ -I-X 3 -I-X 4 -I-X 5 -I-xe), ui = ^{xi-X2), vi = ^{xi + X2 - 2x‘i) 

vS \2 v6 

U2 = -E(a;4 - X5),n2 = ■E(x4X5 - 2 X 6 ), W = ^{xi+X2 + X 5 -X 4 ^-X^-Xq) ( 82 ) 

The transformed Hamiltonian reads: 

TT _ 9"^ , 2.,2 , 2 , 2 , 2 , 2 , 2^ 






^ 9Ai(nf-bnf)2 ^ 9A2(n2 + ^^2)^ , 


9 


+ w‘^ + uf + U 2 + Vi + V 2 2 (uf — 3«ixf)^ 2 (^2 — 3 ^ 2 X 2 )^ 6 w 


2 ■ 


( 83 ) 


This Hamiltonian is not separable in {t,ui,vi,U 2 ,V 2 ,w} variables. We first introduce the following 
coordinates transformation ; 


ui = risin(^i, xi = n cos (/9i, 0 < ri < oo, 0 < (pi < 2 it 

U 2 = r2sin(/92, V 2 = r2Cos(/92, 0 < r 2 < oo, 0 < v ?2 < 27r 

The Schrodinger equation is then written: 


E 

dt'^ 


^2 92 IQ 


dw‘2 5r2 ri dri dr^ r2 8x2 6w‘^ 


d 1 cl g 2/ 9 2 2 2\ 

+ xE + 0J^{ri +rl + e + u;2) 


-F- 


9 




■*■^2 

^2 


■2 -|- rl -I- t 2 -p w'^ rf 


9 A 2 


92 




9Ai 


92 




dip2 2sin2(3(/22) 


d(fi 2sin2(3(/2i) 


'i’{t,w,ri,r2,p>i,(p2) = E'Ii(t,w,ri,r2,p)i,p>2) . 


( 84 ) 


( 85 ) 
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The potential involved in the equation (1851) 

V{t,w,ri,r 2 ,ipi,(p 2 ) = 0 J^{rl + rl+t^+ v?) + 




+ 


g 


1 9Ai 

+ —o „ . o —r + 


has the general form 


rf + rl + ' 6ui2 

l_ 9A2 

rl 2sin^(3(/?i) ' 2sin^(3(/?2) 


/i(<y?i) , 12 (^ 2 ) 


V{t,w,ri,r2,(pi,^2) = f{t,w,ri,r2) H-5-^ 


This suggests the wave function to be factorized as follows 

2 

T(t,u;,ri,r2,(/?i,(/?2) = T(t, u;, ri, r2) n • 

i=i 

The equation (I85p will be solved in two steps. Firstly we solve 

/ £ 


( 86 ) 


(87) 


( 88 ) 


and 


dip"^ 


I ^ 2sin2(3v9i)) 


+ 


9 A 2 


dcp2 2sim(3(^2) 




(89) 


(90) 


on the interval [0,7r/3], with Dirichlet conditions at the boundaries. The Bnj,j = 1,2 are the quan¬ 
tified eigenvalues of the equations ()89I90I) . respectively given by 

= 6W = 3(n, + ^ + a,), n, = 0,1,2,.... j = l,2. (91) 

The associated eigensolutions are given in terms of the Gegenbauer polynomials Cn'^ 

= (sin3(/9j)5+“jCn|"^“'^(cos3(^j), 0 < (/9 < |, =0,1,2,.., aj = I + 2Xj . (92) 


The second step consists in the resolution of the following Schrodinger equation 


d2 q2 q2 


Id 1 d ^ g 


dt'^ dw'^ dr\ ri dri dr^ r 2 dr 2 Qw'^ 


+u;\ri+ri+t^ + w^) + 




rf + r 2 +t‘^ + 


?(1) C?(2) 


I ■’ni I -1^112 
~r o ~r 


2 r (^) ) ^ 2 ) — T^ni,n 2 '^’ni,n 2 


{t,w,ri,r 2 ). 


(93) 
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This Hamiltonian is not separable in {t,w,ri,r 2 } variables. We then introduce the following hyper- 
spherical coordinates : 

t = rcosa, rc = r sin a cos 0, ri = r sin a sin 0 sin/3 r 2 = r sin a sin 0 cos/3 
0 < r < oo, 0<a<7r, O<0<7r, 0</3< 7r/2. (94) 

Since ri,r 2 are positive we have (3 G [0,7r/2]. 


The Schrodinger equation (IU5D becomes : 
\ h d 


+ 


dr'i 

1 

sin^ a 


dr 


+ u;V + 4 + 


52 ^ d 

^ - 4 cot a— 
da^ da 


^ 52 c? g 

QQ2 *201 QQ Q gQg2 Q 


+ 


1 


sin2 9 


52 ^ d 


B: 


( 1 ) 

rtl 


+ 




( 2 ) 

n2 


sin2 /3 cos^ (3 


^ni,n2 ^)/^) — .E'ni,n2^ni,n2 


(r,a,9,^). 


(95) 


This Hamiltonian may be mapped to the problem of one particle in the four dimensional space with 
a non central potential of the form 


V{r,a,e,/3) 


fi{r) + 


1 

9 ■ 2 

sin a 


f2{e) + 



fm 


(96) 


As in the second section, the problem becomes separable in the four variables {r, a, 0, /3}. To find the 
solution we factorize the wave function as follows : 


4- 


ni,n2 


{r,a,6,P) = 




m,ni,n2 


(a) &j,m,ni,n 2 {0)H, 


m,ni,n2 




r'°l‘^ sin2 a sin0^/2 y^sin 2(3 

Accordingly, equation (1951) separates in four decoupled differential equations: 


(97) 


d‘^ {bnif - 1/4 - 1 / 4 ' 

dl3'^ sin2 /3 cos^ (3 

U) 


77m,ni,n2 (/^) — C*m,ni,n2,n2 (/^)) 


in terms of the bii- of Eq. ([9T]) . 


{Cm,ni,n2 4) 9 

d0^ sin2 0 6 cos^ 0 


0j,m,ni,n2 (^) — jj2 0j,m,ni,n2 (^)) 


a ^ ^j,m,ni,n2 4 1 . i _ A f 


and 


, 2 2 , 7^ ^£,?',m,ni,n2 4 

+- 2 -^ 

ar^ 


Bk,£,j,m,ni,n2{^) Bk,i,j,m,ni,n2 d^fc,^,j,m,ni,n2 (^) 


(98) 


(99) 

( 100 ) 


( 101 ) 
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The regular solutions of Ea. (j98p . for /3 G [0,7r/2] (remind that ri,r 2 > 0), with Dirichlet conditions at 
the boundaries read 


H 


m,ni 


= (sin/3)<^+^ (cos/3)^"2+5 

TT 


0 < m = 0 ,l, 2 ,..., 


( 102 ) 

(103) 


in terms of the Jacobi polynomials Pm^’ ■ The associated eigenvalues Cm,ni,n 2 are given by 


C*m,711,712 — ^m,ni,n 2 ’ Cm,ni,n 2 — ( 2 m + + 1 ), m — 0 , 1 , 2 , ... . 


(104) 


Since bn^ / 0,j = 1,2 the Hamiltonian of Ea. dMl) is a self-adjoint operator. 

The second angular equation for the polar angle 9 reads : 

Cm,ni,n2“l/4 - 1/4 

“T . o ^ “T 


d02 


sin^ 9 


cos2 9 


Pj,m,ni,n2 1 ®j,7n,7ii,7i2 (^) — 0) 


(105) 


with Cm,ni,n 2 given by Eg. (11041) and d = < 7/6 -|- 1/4 with the condition that g > —3/2. Since 

,112 / 0) the Hamiltonian of Ea. (|105p is a self-adjoint operator if d / 0 i.e. g / —3/2 


^m,ni 


The regular solutions of Eg. (11051) in ]0, 7 r/ 2 [ with Dirichlet conditions read, in terms of the Jacobi 
polynomials. 




j,m,ni,n2 


(d) = 


0 < 


(sin (cos (cos 29), 

„ TT ,1 I ^ 

9<-, j = 0,l,2,..., d=-^l + ^. 


(106) 

(107) 


The index -|- means that the tc-coordinate is positive. The solutions in the interval ^ < 9 < n are 

0jr4„..„(«) = (sin«)'“'".'”=+i(-cose)"id‘‘”’"‘’""‘'T“2()), (108) 

^ < 9<7r, j = 0,1,2,.... 

The index — means the u^-coordinate to be negative. Note that in the interval ] 7 r/ 2 , 7 r[ we have 
— cos d > 0 and sin d > 0 so that the real power of these positive values are defined. 


In both cases, the eigenvalues Dj^rn,ni,n 2 of Ea. (ll05p are given by 


D = d^ 

^j,m,ni,n2 '^j,m,ni,7i2 > 


d. 


'j,m,711,712 


— (2j “1“ Cm,m,712 T d -|- 1), j — 0,1, 2,.... 


(109) 


The regular eigensolutions and corresponding eigenvalues in the interval ] 0 , 7 r[ of the third equation 
poop read, respectively. 
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Ge,j,m,ni,n 2 {a) = (sin '"2 + 2 ^ ^ (cOS Of) , ^ = 0,1,2,..., 

2 f 1\ 

^£,j,m,ni,n2 ~ ^£,j,m,ni,n2^ ^£,j,rn,ni,n2 ~ {^~^dj,m,ni,n2 2 ) ’ -^ = 0 , 1 , 2 ,. 


Note that the choice dj^rn,ni,n 2 > 0 implies for every value of {j, m, ni,n 2 } the function 
to vanish at the boundaries of the interval ]0,7r[. 


Finally, the reduced radial equation reads 


dr"^ 


+ + 




_ 1 
J,m,711,712 4 


d^k,£,j,m,ni,n2 1 J,m,ni,n2 (^) 0. 


Taking into account the definition of ^£j,m,ni,n 2 J Ea. (|llip . the constraint 

/ . 11\2 
fi + ^rj,m,ni,n 2 ~ T ( ^ T + 2m + 3ni + 3n2 + d + 3oi + 3o2 H —— \ >0 

Vni > 0, Vn 2 > 0, Vm > 0, Vj >0, > 0, 

happens, for every positive value of /r. Since /r + vlrj,m,n is minimal for re = 0, j = 0, m = 0, ^ = 
ai = 02 = d = 0, it put constraint on negative values of /r, which has to satisfy 


'VTT^ 


< /i < 0 


We introduce the auxiliary parameter 


j,m,ni,n2 


defined by 


f^£,j,m,ni,n 2 — + ^rj,m,ni,n 2 — y + 2m + 3rei + 3re2 + d + 3ai + 3a2 + —^ 

The solution of the radial equation (111211 vanishing at r = 0 and r ^ oo are : 


Fk,£,j,m,nun2{r) = 1 - 2 +^ exp 4"^—i-2)(a;r2), 


k = 0,1,2. 


in terms of the generalized Laguerre polynomials and the associated eigenenergies are given by 
Fk,£,j,m,ni,n2 — 2oj{2k T K£J^m,ni,n2 T f); ^ — 0) fj ^ — 


( 110 ) 

( 111 ) 

,n2 (®) 

( 112 ) 

(113) 
0 and 

(114) 

(115) 

(116) 

(117) 


19 









Collecting all pieces, we conclude the physically acceptable solutions of the Schrodinger equation ([8^ 
to be given by 


X (sin Q,)2i+2m+3ni+3n2+d+3ai+3a2 + |(^j2i+2m+3ni+3n2+d+3ai+3a2 + ^) 

X (sin 0)2m+3ni+3n2+3ai+3a2+3^gg l pj2m+3ni+3n2+3ai+3a2+4.d) ^0) 

x(^)(sin/3)3’^l+3ai + f^^Qg^)3n2+3a2 + fpj^3«l+3“l + i3’"2+3a2+f)^^^^ 

V2 

xs( 7 n(sin | sin (cos 3<y9i), 

xsgn{sin | sin 3(/92r^'''2 ^ (cos 3(p2), 

with e = 1 for bosons and e = — 1 for fermions and 

k = 0, 1,2, ...,i = 0,1,2, ...,j = 0,1,2, ..,m = 0,1,2, ...,ni = 0,1,2, ...,n 2 = 0,1,2, 

vr ^ vrl — esvr „ 3 — e^vr 

0 < ,O</3<-,^-^-<0<^-^-,e3 = ±l,O<a<7r 


(118) 


1 


ai — —\/l + 2Ai 


1 


0-2 — — \/l + 2 A 2 , d — — \/1 H—— 


29 


It has to be noticed that, for Bose statistics, a 6 pathology occurs in (|118l) for aj = 1/2 {Xj = 0,j = 1, 2) 
and d = 1/2 (g = 0) 

The normalization constants -/Vfc,rj,m,ni,n 2 can be calculated from 


r+oo 


r^dr / sin^ a da 


r(3-e3)7r/4 
/ ( 1 - 63 ) 71/4 


7r/2 — — 

sin3 0 dO [ sin(2/3) d/3 d(/?i dip 2 


^ k,i,j,m,nin2 (^) ®, d, /3, (/>! , ^^ 2 )^ k',j' ,m' ,n'-^,n2 iTi 1 1 ¥^2) — ^k,k'^l,l'^ni,n[ ^n2,n'2^k,i,j,m,n\,n2 ■ 

As above, use is made of the orthogonality properties of Gegenbauer, Jacobi and Laguerre polynomials. 

The eigenenergies are 




711^712 


— 2cj I 1 -\- ‘2h \ jjj -\- 2ij ‘2 tti -\- 3/1]^ -t- 3712 d 3cii H- 3ci2 H —\~~ 


11 


(119) 


Partial degeneracies are observed i.e. all solutions for which i + 2j + 2m + 3ni + 3n2 = Ai is satished, 
N, being hxed. 
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5 Conclusions 


In this paper, we have proposed and exactly solved few-body quantum problems of four, five and six 
particles in the D = 1 dimensional space. The particles are confined in a mean field generated by 
harmonic or attractive Coulomb-type forces. The interactions between particles are governed by two- 
body inverse square Calogero potentials with both translationally and non-translationally invariant 
many-body forces, for the three problems treated. The difficulties encountered concerning the explicit 
expression of the eigenfunctions, in problems concerned by the pairwise Calogero interaction between 
each pair of particles, have been overcome by restricting the number of two-body interactions. Our 
choice singularizes the interaction between one cluster of three particles with one particle for the four- 
body problem. For the five-body problem , we singularize the interaction between two clusters of two 
and three particles respectively. Finally, for the six-body problem, we choose the interaction between 
two clusters each one formed by three particles. 

Appropriate coordinates transformations, lead to the solution of the stationary Schodinger equation, by 
providing the full regular eigensolutions with the eigenenergies spectrum for these few-body bounded 
systems. 

For the four-body problem we pay attention to the irregular solutions. The coupling constants domains 
for which the irregular solutions, become square integrable are determined. Also for the four-body 
problem we have replaced the harmonic confinement by an attractive ” Coulomb-type” interaction, 
and the solutions are given for the bound states only. This study can be extended straightforwardly 
to the five and six-body problems. 

The present results suggest to extend the construction of exactly solvable models to larger systems of 
particles. 
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